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Abstract 

The time development of the price of a financial asset is considered by constructing and 
solving Langevin equations for a homogeneously saturated model, and for comparison, for a 
standard model and for a logistic model. The homogeneously saturated model uses coupled 
rate equations for the money supply and for the price of the asset, similar to the coupled rate 
equations for population inversion and power density in a simple model of a homogeneously 
broadened laser. 

Predictions of the models are compared for random numbers drawn from a Student's 
t— distribution. It is known that daily returns of the DJIA and S&P 500 indices are fat 
tailed and are described well by Student's t— distributions over the range of observed values. 
The homogeneously saturated model shows returns that are consistent with daily returns for 
the indices (in the range of —30% to +30%) whereas the standard model and the logistic 
model show returns that are far from consistent with observed daily returns for the indices. 

Keywords: asset prices; returns; homogeneously saturated; logistic; standard model; 
Langevin; Student's t— distribution; truncation 
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i 1 Introduction 



2 The standard model for the development in time of the price of an asset is geometric Brownian motion 

3 with a deterministic growth (or drift) rate [U [2j [3] , jU Ch 16.4]. The region of support for the normally 

4 distributed noise driving term of the Brownian motion is taken as — oo to +oo, which allows for infinite 

5 prices. This choice of support does not cause difficulty in the standard model as the probability of a large 
e price is essentially zero; the probability of a noise driving term with value x goes as exp(— x 2 ). However, 
7 the choice of a 'standard' model, which includes an infinite region of support for the noise driving term, 
s makes pricing assets and pricing options based on these assets difficult when the underlying probability 

9 density functions (pdf ) is a fat tailed distribution rather than a normal distribution. 

10 It is known that daily returns of the DJIA and the S&P 500 indices are described by Student's 

11 t— distributions [31 [5] and that Student's t— distributions have fat tails. Hence prices predicted by the 

12 'standard' model with a fat tailed distribution for the noise driving term can be very large (essentially 

13 infinite) with a non-zero probability. Here the standard model is defined to include the same region of 

14 support of — oo to +oo for the noise driving term as for Brownian motion, and the returns are (naively) 

15 assumed to follow the distribution that fits the observations over the the full region of support. For 

16 the standard model with an underlying fat tailed distribution, integrals to price European call options 

17 diverge [6J U\ E] • Thus a standard model with a fat tailed distribution is not adequate, and a standard 
is model with a normal distribution does not match observed daily returns. 

19 There exist several approaches to price options and assets when the underlying distribution is a fat- 

20 tailed distribution. One approach is to modify the tails of the distribution such that the contributions far 

21 into the tails are negligible while not affecting significantly the central portion of the distribution, which 

22 fits well the observed data. This can be accomplished by multiplying the distribution by an exp — \a\t 2 

23 envelope function [8] or by using a generalized t— distribution that has terms of exp — |a|t 4 [HI [10]. It has 

24 been demonstrated that an exp — \a\t 2 envelope multiplication is obtained for a Student's t— distribution 

25 by truncation of the volatility in the chi-normal mixture that leads to a Student's t— distribution [TT1IT2"] . 

26 It is possible to price assets and options with a standard model that uses fat tailed distributions by 

27 capping the value of the asset or by truncating the underlying pdf. Using the arbitrage theorem and the 

28 constraint that the process be fair (i.e., the price is a martingale), prices for European call options can 
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29 be found for fat-tailed distributions if the distribution is truncated or if the value of the asset is capped 

30 [13]. Truncation or capping keeps integrals, which are needed to price options, finite [3 [TJ [5]. 

31 Truncation of the pdf for the returns and capping the value of an asset might seem to be ad hoc 

32 solutions that allow for any price [HI [7]. Truncation is based on conditional probability, which is a sound 

33 mathematical concept. In truncating, one accepts that there is a probability that the value of the asset 

34 will exceed the truncation. The option writer has the ability to select the risk that the value of the 

35 underlying asset will exceed the truncation and to price the option for the selected risk. 

36 For capping, both the writer and buyer must agree that the price of the option is based on a maximum 

37 value for the asset. Both sides should recognize that there is a finite probability that the value of the 

38 asset will exceed the agreed upon maximum value but accept that this maximum value will be used. 

39 This controls the risk, keeps the integrals finite, and allows for the options to be priced. 

40 A different approach to price assets and options when the underlying distribution is a fat-tailed 

41 distribution is to allow for saturation of the price of an asset by depletion of the resource that supports 

42 the price (i.e., by depletion of the reservoir of money that is available to purchase the asset). This is 

43 the approach that is investigated in this paper. In this paper a homogeneous saturation model for the 

44 price of an asset is constructed and compared to the standard model and a logistic model to gain insight 

45 into the pricing of financial assets and options based on these assets. Random numbers drawn from a 

46 (fat-tailed) Student's t— distribution are used to compare the predictions of the models. 

47 The homogeneous saturation model borrows from laser physics, wherein coupled rate equations are 

48 used to describe the interaction between the output (equivalent to price in the pricing of assets considered 

49 here) and inversion (equivalent to the reservoir of money available to purchase the asset). The saturation 

50 in a laser keeps the output finite for finite input and ensures that the power output equals the power 

51 input. The saturation in the pricing of an asset keeps the price finite for a finite supply of money, such 

52 that the integrals that are required to price an option based on the asset remain finite. 

53 A homogeneously saturated model for the development in time of the price of financial assets is 

54 investigated in this paper, and is compared to a standard model and a logistic model. Simple Langevin 

55 equations for the time development of the price of an asset, which are first order differential equations 

56 with noise driving terms and which should be interpreted as integral equations [HI pg 172] (U Ch 10.2], 
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57 are constructed and solved for the standard model, a logistic model, and for a homogeneously saturated 

58 model. The predictions of the models are compared for fat tailed noise driving terms. The logistic 

59 model is the standard model with a non-linear saturation term whereas the homogeneous saturation 
eo model follows laser physics [151 ESJ 02] • I n a logistic model, the non-linear saturation term keeps, e.g., 
ei the voltage finite in a Van der Pol oscillator [T5l Eq 6, pg 46] and populations finite in competitive 

62 environments JTSJ Q3J]. For the homogeneous saturation model, coupled rate equations for the price of 

63 an asset and the money supply supporting the asset are postulated and solved. This coupling between 

64 the price and the money supply keeps the price finite and allows for pricing of an European option with 

65 fat tailed distributions. Similarly, the coupling between the inversion and output in a laser keeps the 
ee output of the laser finite and equal to the input. Since only one asset is considered in the simple approach 
67 presented in this paper, the money supply is saturated uniformly (homogeneously) by the asset. It is 
ea possible to envision multiple assets interacting with the money supply and with other reservoirs. In this 

69 approach one might allow for inhomogeneous saturation, or homogeneous saturation, or some mixture 

70 of the two limiting cases of saturation, and allow for low prices to stimulate purchases, large prices to 

71 stimulate sales, and spontaneous decisions. 

72 Toth et al. [20] and Bouchaud et al. [21] studied the order book and described the microstructure 

73 of the market. These researchers reported that the stability of markets depends on a precise balance 

74 between supply and demand, and that price is a steady state and not an equilibrium. Smith et al. [22] 

75 used a rate equation for the density of the order book to understand how prices depend on the rate of 

76 flow of orders. The homogeneously saturated model presented in this paper is consistent with the ideas 

77 of these researchers. The homogeneously saturated model is a phenomenological approach, as is the 

78 work of these authors, that uses coupled rate equations to find a balance between supply and price. The 

79 emphasis in this paper is on simple models for the price of a financial asset and not on understanding 
so the microstructure of the market. 

si Bouchaud and Cont [23] and Bouchaud [24] developed a phenomenological Langevin approach to 

82 study market crashes. They developed a second order differential equation (DE) that is linear in de- 

83 mand minus supply but non-linear in price and solved this DE for the time rate of change of the price 

84 of the asset. The DE for the time rate of change of the price contained a logistic type saturation term. 
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85 The derivative of the price went to — oo once the price of the asset exceeded a threshold value. It was 

86 concluded that crashes are the result of a succession of improbable and unfavourable events, that no 

87 precursor to market crashes exists, and that the market behaves as an adaptive system. The homo- 
ss geneously saturated model presented in this paper is a phenomenological model and is composed of 

89 first order coupled differential equations. The feedback inherent in the coupled rate equations makes 

90 the output of the homogeneously saturated model adapt to changes in price and money supply. The 

91 emphasis of this paper is on three simple models for the development in time of the price of an asset 

92 and not on market crashes. 

93 Grassia [25] added market delay and feedback to obtain a linear second order DE for the price. 

94 Grassia also solved for the time rate of change of the price and thus obtained a Langevin equation. 

95 Grassia studied the time dependence and stability of the price. He found that quenching ensured long 

96 term bounding of the price of the asset. Richmond and Sabatelli [26] developed a Langevin model of 

97 interacting agents to understand fluctuations of the prices of financial assets. They used their results to 

98 understand the personal incomes of several countries. 

99 Anteneodo and Riera [27] developed a non-linear mean reverting Langevin model to study the 

100 stochastic dynamics of volatilities. Anteneodo and Riera showed that additive-mutliplicative processes 

101 are required to obtain fat tailed distributions. In this work the underlying fat tailed distribution is 

102 accepted as a fact [281 [291 E EH] an d is used with the three models to investigate the pricing of assets. 



104 In the equations that follow, S(t) is the value of an asset at time t, S = S(0) is the value of the asset 

105 at t = 0, a is a drift rate, a is a scale parameter, and f(t) is a stochastic process. 
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2 Standard Model 
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The standard model for the time development of the value of an asset is 




a S{t) + a S{t) fit) 



(1) 
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with solution 
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-t 

S(t) = S exp / (a + af (77)) cfy . (2) 



R(t) = Ini^P-) = [\a + af( V )) d V . (3) 



108 The return R(t) is 



S 

109 If the integral of the noise wanders to infinity, which is possible with fat-tailed distributions, then 
no the return goes to infinity. The equation for the return shows that the noise contributes to the return 
in and dominates when a f(t) > a, which happens routinely since a is typically a small number. 
112 The equation for the time development of the average value for the asset is given by 

d S(t) = (a+^) S(t) (4) 



dt w V 2 

with solution 



S(t) = S e( a+ ^ t . (5) 

114 The average value tends to infinity as t tends to infinity for 2 a + a 2 > . 

115 The equation for the time development of the average value of the square of the value of the asset is 

±S2(t) = 2 (a + a 2 ) S*{t) . (6) 

lie These equations give the variance of the value of the asset given that the asset was worth S Q at t — 
117 as 



Var(S(t)) = S 2 e ( 2a + <j2 )* (> 2 * - l) , (7) 

us which is the variance for a log-normal distribution. This is expected, as the equations for the average 

119 values of S(t) and S 2 (t) are obtained using a Langevin approach with (f(t)) = 0, (/(ii)/(^2)) = — ^2), 

120 and higher order expectations = 0. The same results for the average values are obtained by use of Ito's 
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121 calculus |U pg 189]. 

122 As t approaches infinity, the variance approaches infinity if a + a 2 > 0. 

123 It has been assumed that the conditions for the Langevin approach (i.e., that the diffusion coefficients 

124 Dk = for k > 2) hold [U Ch 10] for equations that are driven with noise sources that are distributed 

125 as Student's t— distributions, or, equivalently, that Ito's calculus holds. It is known that Student's 

126 t— distributions fit market returns and that the time correlations for daily returns are approximately 

127 delta function correlated. Thus the assumption is well motivated and seems reasonable in that it is 

128 consistent with observations. 

129 In the following sections the differential equations for the time development of the value of an asset 

130 for the logistic and the homogeneous saturation models are solved. The solutions are then used with 

131 Student's ^-distributions to investigate prices of assets. 



132 3 Logistic Model 

133 The time development of the value for an asset using a logistic model is 

±S(t) = aS(t) - (3S 2 (t) + aS(t)f(t) (8) 

134 with solution 

C p Jo a+ (7/(77) dT? 

S(t) = ° . c . (9) 

1 + f3S J"* etf" + a* dC 

135 Note that the standard model is obtained from this logistic model in the limit (3 = 0. 

136 The equation for the development in time of the average value of an asset that follows the logistic 

137 equation above is 

±s{t) = (« + y) s(t) - 0s\t) (io) 

138 with solution 
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S(t) 



S (a + 



ps 



a + 



- 1 



139 The average value of the asset remains finite for all time and approaches the value 



lim S{t) 

t— >oo 



a + 



P 



(12) 



wo as t tends to oo. 

141 Figure 1 is comprised of plots of the average value of S(t) as a function of time. The limiting 

142 behaviour as t approaches infinity is clear in the figures. For finite a + a 2 /2 the value of S(t) saturates 

143 to a finite value for large t. This behaviour is distinct from a random walk where J a + af(t)dt wanders 

144 or drifts to infinity. The standard model is a random walk and is obtained from this logistic model by 

145 setting (3 = 0. 
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Figure 1: Time development of a logistic variable with finite drift for four different values of the saturation 
parameter. 



The equation for the time development of the square of the value of the asset is 

J 3 

^-S\t) = 2(a + cr 2 )S 2 {t) - 2f3~S~ 2 \t) 



(13) 



with solution 
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S\t) = S " ^ + ° r 2 • (14) 

(e-<" + » , > 1 (^S„ - a - <7 2 ) - 3S„) 2 

148 In the limit as t approaches infinity, the variance of S(t) under a logistic model remains finite and 

149 equals 

a + a 2 + 3 (eff 
lim Var(5(t)) = V 7 . (15) 

150 An approximation for Eq. (9) can be obtained. The outer integral in the denominator of Eq. (9) 

151 presents a challenge, as the value of the integral depends on the values that the stochastic process takes 

152 for each point in time. If one models the stochastic process as small steps in the same direction and 

153 of the same magnitude, then the integral can be converted to a summation and this summation can be 

154 evaluated. In this simplified model, which is adopted to allow quick evaluation of the value, the picture 

155 then is one of noise as small steps in the same direction, not as abrupt, large jumps. 

156 Let 

af( V ) dr, = W{t) (16) 

o 

157 and let the Weiner process W(t) = W. Under this simplified picture of the stochastic process, the 

158 integral in the denominator of Eq. (9) can be approximated as 

+ d C «± J2e^> = il^r. (17) 

</o iV i=0 iV 1 — e n 

159 Provided that the interval [0, t] is subdivided into N intervals such that N » at + W, then a 
wo series expansion of the exponential in the denominator can be used to find that 

rt p<Xt + W _ 1 







a + W/t 



and 



So (« + f ) e at + w 
S{t) ~ a + f + (3S {e«t + w _ I) • ( 19 ) 
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162 For large at + W, 



lim S(t) = (20) 

at + W^-oo (3 

163 and the value of the option approaches infinity as a + W/t approaches infinity. However, the approach 

164 is not exponential, as is the approach for the standard model, Eq. (2). The return is not linear in a + 

165 W/t; the return is the logarithm of a + W/t. W/t is the average step size (i.e., the total change owing 
lee to noise divided by the time taken to make the change). 

167 Figure 2 is a plot of the S(x) versus x for x = a + VF(1). S(x) defined in this manner is the value of 

lee the asset for one day later (i.e., for t — 1) with x the value for the accumulated drift and noise over the 

169 one day. The initial value of the asset, S , was taken as 50. — gives the standard model. The three 

no non-zero values of /3 of 0. 05/50, 0.01/ S , and 0. 02/ S Q are the same as in Table 1. The logistic model 

171 provides some saturation for large values of the drift parameter a and the accumulated noise W(l) and 

172 no saturation for negative values of the accumulated noise. It is interesting to contrast the behaviour 

173 as a function of x = a + W(l) with the time development of the solution to a logistic equation. For the 

174 time development of the solution to a logistic equation, as t approaches infinity, the value of the asset 

175 approaches x/(3. For a given x that is constant in time, the logistic equation saturates. However, if as 

176 shown in Fig. 2 x is not limited then the solution to the logistic equation is not limited. 

8000 
6000 

S4000 

2000 

-20 -10 10 20 

x 

Figure 2: Value after one day as a function of the drift plus noise accumulated over one day for different 
values of the saturation parameter. A logistic equation was solved to find the value S(x) 

177 Figure 3 shows the return as a function of x where the return is calculated as the natural logarithm 
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178 of S(x)/S . Since S(x) was calculated for t =1, the return is the daily return given the integral of the 

179 drift and noise over a time frame of one day. The logistic equation shows some saturation of the positive 

180 return but no saturation for a loss. 

8 

4 

3 

a> 

-4 



-20 -10 10 20 

x 

Figure 3: Daily return predicted by a logistic equation for various values of the saturation parameter 
and a function of the drift and the noise accumulated over one day. 

181 Table 1 lists descriptive statistics for simulated values. Each cell of results contains the ordered pair 

182 S(l), R{1) where S(l) is the value of the asset after one day and R{1) is the daily return, defined in 

183 Eq. (3) as \n(S (1) / S (0)) = \n(S (1) / S Q ) . To illustrate the effect of the saturation, samples were drawn 

184 from a Student's ^-distribution with v — 2 degrees of freedom multiplied by 10xcr/sqrt(365) where a is 

185 an annualized volatility of 0.3. Fits to the daily returns for the DJIA and the S&P 500 show that v — 3 

186 would be appropriate [5]. A drift parameter a = 0.15/365 was assumed. The initial value for S(t) was 

187 taken as 5" = 50. 

188 Descriptive statistics for 4096 samples of W(t) drawn from a v = 2 Student's t- distribution and 

189 scaled by 10xcr/sqrt(365) are maximum value = 7.70, minimum value = —6.01, mean = 0.007, standard 

190 deviation = 0.43, and kurtosis = 86. 

191 Table 1 Descriptive statistics for simulated values after one day, S(l), and daily returns, 

192 using Eq. (19) for 4096 samples drawn from a Student's ^-distribution with v — 2 

193 degrees of freedom, a = 0.0041, a = 0.157, and S a = 50.0 




11 



parameter 


p = o 


PS Q = 0.05 


ps Q = 0.1 


PS = 0.2 


max 


lxlO 5 , 7.70 


7210, 4.97 


3726, 4.31 


1895, 3.63 


min 


0.12, -6.00 


0.12, -6.01 


0.12, -6.02 


0.12, -6.03 


mean 


104.3, 0.008 


55.9, -0.04 


51.1,-0.092 


45.4, -0.18 


std dev 


2139, 0.43 


170, 0.39 


93.4, 0.38 


50.1, 0.36 


kurtosis 


2237, 86 


1386, 55 


1121, 53 


879, 53 



195 The logistic equation shows some limiting of the range of returns for positive returns. Note that the 

we minimum value of the daily return, R(l), is the same as W(t) (cf. the minimum returns for non-zero 

197 P to the minimum return for P — and to the minimum value for the descriptive statistics for 4096 

198 samples, found in the paragraph preceding Table 1). The minimum return seems to be unaffected by 

199 the logistic equation. 

200 For the standard model (P = 0) the values of an asset for one day later, S(l), lie in the interval from 

201 0.12 to 1 x 10 5 with daily returns, R(l), in the interval of —600% to 770 %. For a logistic model for the 

202 time development of the value of an asset, S(l) lies in the interval of 0.12 to 3276 with daily returns 

203 of —602% to 431 %. Returns of < —100% make little sense. The logistic model reduces the maximum 

204 value of the asset owing to noise. The logistic model also affects the mean value. Note that the mean 

205 reduces from 0.008 to —0.18 as P is increased from zero. 

206 4 Homogeneous Saturation 

207 The logistic equation can be justified from a rate equation analysis. 

208 Let M(t) be the amount of money that is available to invest in an asset. Let N be the rate at 

209 which money is pumped into the reservoir of money M(t) that can be used to purchase the asset and 

210 let PxM(t)xS(t) be the rate that money is removed from the reservoir owing to purchases of the asset. 
2n Let r be a characteristic time constant that allows for money to be removed or added to the reservoir, 

212 depending on whether M(t) is greater than or less than some value M Q . At this level of discussion the 

213 value of M Q is immaterial as M /t can be combined with N . 

214 A rate equation for M(t) is then 
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4-M(t) = N — f3 x S(t) x M(t) + ^ — . (21) 

at r 

215 All parameters in the rate equation have a time dependence, but it is assumed that these parameters 

216 change slowly in time owing to inertia in the system and to gradual evolution of tastes and performance 

217 with time. Thus each point in time is assumed to evolve about a steady state but it is accepted that 

218 this steady state point will also evolve slowly. In the concepts of Lax [U pg 162], fluctuations about 

219 an operating point are considered, but the operating point is a point of steady state and not a point of 

220 equilibrium. In steady state, the time derivative equals zero and 

- f^A) = TTJsm • (22) 

221 The second form for M(t) has been recast so that leading terms in a series expansion of 



M(t)xS(t) = aS(t) - (3S 2 (t) + — S 3 (t) - ^ SHt) + ... (23) 



222 give a logistic equation for the time development of an asset when 

±S(t) = M(t)S(t) + aS(t)f(t) . (24) 
at 

223 A closed form solution for the equation above for S(t) using the full form for M(t), Eq. (22), does not 

224 appear to exist. 

225 A noise term is traditionally added to the value of the stock. However, noise can be added to the 

226 rate equation for the money available to invest in an asset. This approach seems fundamental in that 

227 it is the amount of money chasing an asset that determines the value of an asset. The value of a liquid 

228 asset is a visible and easily obtained attribute but perhaps not a fundamental quantity. 

229 If one adds a noise term to the rate equation for M(t), Eq. (21), then in steady state 

M(t) - + */(') _ (" + */(*)) (25) 
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230 In the last form for M(t) the symbols have been redefined. The equation for the time development 

231 of the value of an asset becomes 



±S{t) = M(t)S(t) 



aS(t) + aS(t)f(t) 

1 + ps(t) 



(26) 



232 



In this approach the noise is ascribed to fluctuations in the amount of money available to invest in 



233 the asset. The value of the asset under this picture is 



S(t) 



S e at + w ^ 



(27) 



e /3(S(t) - So) 



e P(S(t) - So) 



234 a transcendental equation that can be solved for S(t) given a, (3, S Q , and W(t). The simple dependence 

235 on W(t) makes simulation straight forward. 



237 of 0. The case /3 — corresponds to the standard model, Eq. (1). The saturation denominator limits 

238 the range of values for the asset S(t) and for the return R(t). Note that for the case f3 = 0, S(20)/So 

239 = e 20 = 5xl0 8 and 5(-20)/5 = e" 20 = 2xl0~ 9 , whereas for the case (3 = 1.0, S(20)/S = e a33 = 

240 1.39 and S(— 20)/ Sq = e -0,49 = 0.61 Clearly the homogeneous saturation effectively limits the return: 

241 the standard model, /3 = 0, predicts increases of 500 xlO 6 times the initial value of the asset whereas 

242 a homogeneous saturation model with (3 = 1 predicts an increase of 1.39 times the initial value of the 

243 asset, both for an accumulated drift and noise of x = 20. 

244 The descriptions for M(t) and S(t) shown in Eqs. (25) and (27) are similar to a simple description 

245 of a homogeneously broadened laser amplifier [17] . It is possible to find analytic (albeit transcendental) 

246 solutions for the power in each mode in the case of a simple description of a homogeneously broadened 

247 gain medium [T7]. Thus, it might be possible to solve for multiple assets interacting with the same 

248 reservoir of money for different a and different /3 for each asset. 

249 Table 2 lists descriptive statistics for simulated values. Each cell of results contains the ordered pair 

250 S(l), R{1) where S(l) is the value of the asset after one day and R{1) is the daily return, defined in 

251 Eq. (3) as hi(S (1) / S (0)) = ln(S(l)/S ). To illustrate the effect of the saturation, samples were drawn 

252 (as for the results presented in Table 1) from a Student's ^-distribution with v = 2 degrees of freedom 
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Figures 4 and 5 show S(x) and R(x) for as a function of x = at + W(t) with t = 1 for four values 
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Figure 4: Value after one day, S(x), as a function of the drift and noise accumulated over one day, x, 
and for various values of the saturation parameter, 0. (3 = corresponds to the standard model. The 
value S(x) was determined by solving a homogeneous saturation equation. 
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253 multiplied by 10x<r/sqrt(365) where a is an annualized volatility of 0.3. Fits to the daily returns from 

254 the DJIA and from the S&P 500 show that v = 3 would be appropriate [5]. A drift parameter a = 

255 0.15/365 was assumed. The initial value for S(t) was taken as S Q = 50. 

256 Table 2 Descriptive statistics for simulated values after one day, S(l), and daily returns, 

257 R(l), using Eq. (27) for 4096 samples drawn from a Student's ^-distribution with v = 2 

258 degrees of freedom, a = 0.0041, a = 0.157, and S Q = 50.0 



parameter 


= 


= 0.25 


= 0.5 


= 1.0 


max 


lxlO 5 , 7.70 


79.0, 0.46 


64.9, 0.26 


57.6, 0.14 


min 


0.12, -6.00 


28.3, -0.57 


38.5, -0.26 


44.1, -0.13 


mean 


104.3, 0.008 


50.0, -0.0001 


50.0,-0.0002 


50.0, -0.0001 


std dev 


2139, 0.43 


1.58, 0.031 


0.82, 0.016 


0.42, 0.008 


kurtosis 


2237, 86 


88, 79 


86, 76 


86, 79 



260 The data in the Table 2 clearly shows that the saturation included in Eq (27) limits the effect of 

261 the noise on the value of S(t). With the standard model (0 = 0) the maximum value of S(l) is lxlO 5 

262 with a maximum daily return of 770% whereas with = 0.25 the maximum value of S(l) is 79.0 with 

263 a maximum daily return of 46%. Thus saturation tends to limit the range of daily returns. 

264 It is interesting to note that over the history of the DJIA and the S&P 500 the daily returns are in 

265 an approximately symmetric interval of —0.30 to +0.30 [5]. The entries in Table 2 for non-zero show 

266 a truncation of the returns in an approximately symmetric interval. Given the maximum magnitudes 

267 of the returns for the DJIA and the S&P 500, it appears that a value of = 0.5 would be appropriate. 

268 It is also interesting to note that the homogeneous saturation does not affect the mean value, unlike 

269 the logistic saturation that does decrease the mean value as the magnitude of the saturation parameter 

270 is increased. 

271 An approximation to the expression for S(t) can be made by substituting S(t) = S Q exp(at + W(t)) 

272 and making a Taylor series approximation for the exponential term in the denominator. The result is 



S(t) 



p /3(5(t) - So) 



S e at + w ^ 



1 + 0S O e at + w ® - 0S O 



(2i 
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273 In this approximation, for large at + W(t) the value of S(t) approaches 1//3. The approximation for 

274 S(t) above satisfies the differential equation (DE) 

^-S(t) = (a + a x /(*)) x S(t) x (1 - /3 x S(t)) . (29) 

275 In this case the saturation term saturates both the noise and the drift term. The form of the saturation 

276 suggests the first term of a series expansion of a saturation denominator. The solution to this DE has 

277 interesting behaviour. The standard model is obtained for (3 = 0, whereas for (3S Q = 1, S(t) = S(0) = 

278 S Q for all t, and for (3S a > 1 and at + W(t) < 0, negative values of S(t) are obtained. Thus < (3S Q 

279 < 1. 

2bo Table 3 gives descriptive statistics for simulations. Each cell of results contains the ordered pair 

281 S'(l), -R(l). Similar parameters were used in the simulations for Table 3 as were used in the simulations 

282 for Tables 1 and 2. It is interesting to note that the approximate formula, Eq. (28), does not yield 

283 a symmetric interval for the daily return. The data shown in Table 3 show that the minimum return 

284 is not greatly affected by the saturation term whereas the maximum return is greatly affected. Thus 

285 whereas Eq. (28) might be simple and computationally efficient, the equation might not be an adequate 

286 description of losses. 

287 Table 3 Descriptive statistics for simulated values after one day, S(l), and daily returns, 

288 i?(l) , using Eq. (28) for 4096 samples drawn from a Student's ^-distribution with v = 2 

289 degrees of freedom, a = 0.0041, a = 0.157, and S = 50.0. 



parameter 


p = o 


13 S = 0.4 


f3S = 0.8 


f3S = 0.9 


max 


lxlO 5 , 7.70 


125, 0.92 


62.5, 0.22 


55.6, 0.11 


min 


0.12, -6.00 


0.21, -5.49 


0.61, -4.41 


1.2, -3.72 


mean 


104.25, 0.008 


50.5, -0.001 


49.73,-0.011 


49.77, -0.007 


std dev 


2139, 0.43 


9.6, 0.25 


3.8, 0.14 


2.6, 0.10 


kurtosis 


2237, 86 


125, 117 


50, 400 


138, 658 
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291 5 Discussion 



292 The standard model for the development in time of the value of an asset is S(i) = A(t) expcx j£(t)dt 

293 where £(t) is a random variable, a is a positive constant, and A(t) includes the drift. This model predicts 

294 an infinite value for the value of an asset if the accumulated value of the noise a J £(t)dt approaches 

295 infinity. Clearly an infinite value for an asset is not physical [61 [7]. The amount of wealth is limited. 

296 Thus the standard model is valid only for small f£(t)dt and the standard model should not used to 

297 predict values of assets for large J £(t)dt . 

298 Unfortunately, daily returns are known to have fat tails; Student's ^-distributions are known to fit 

299 daily returns well [21 pg 88], [SJ [12] • In contrast to a normal pdf, the probability of a large return is 

300 non-zero for a fat tailed pdf such as a Student's t— distribution. 

301 To price an European call with the standard model and assuming that £ = J £(t)dt follows a Student's 

302 ^-distribution with 3 degrees of freedom, it is necessary to evaluate an integral of the form 

eXp((T ° X 3 5 (30) 

' OH) 

303 where is the value of the strike at time T and At is the expected value of the asset at time T [5]. 

304 The integral equals infinity for a > 0. The exponential numerator from the standard model of the price 

305 of an asset dominates the power law of the fat tailed distribution. Thus an infinite value for a European 

306 call option is found with fat tailed distributions and the standard model [6j [7J El [13]. Figure 6 shows on 

307 a logarithmic scale the value of the numerator as dotted lines and the value of the integrand for < £ 

308 <100. Clearly as £ tends to infinity, the value of the integrand tends to infinity. 

309 A logistic and a homogeneous saturation model were considered as descriptions for the value of an 

310 asset. Simulations show that the logistic model limits, as compared to the standard model, the positive 
3n return whereas the homogeneous saturation model limits both positive and negative returns. Both 

312 positive and negative returns are observed to be limited (i.e., to fall within a range) for the DJIA and 

313 for the S&P 500 indices |5] see, e.g., Fig.l ] 

314 The nonlinear saturation provided by the logistic and by particularly the homogeneously saturated 

315 model limit the value of an asset such that the integral required to price an asset remains finite. The value 
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Figure 6: Plot of the numerator (broken line) of the integrand and the full integrand (solid line) as 
a function of the variable of integration for an integral required to price an European call option as- 
suming the return follows a standard model with a fat-tailed noise driving term given by a Student's 
t— distribution with three degrees of freedom and a standard deviation of \/3cr. The thin tic marks give 
the one-tail critical values for probabilities of 0.99, 0.999, 0.9999, and 0.99999 for the t— distribution. 

316 of the asset S(x) appears to increase linearly with x for the homogeneously saturated model as opposed 

317 to the exponential increase for S(x) for the standard model (c.f. Fig. 4). Thus the integral required 

318 to price an option, Eq. (30), remains finite for the value of the asset as predicted by a homogeneously 

319 saturated model. 

320 A price for an European call option can also be found with the standard model if the value of the 

321 asset is capped or if the distribution is truncated. The thin tick marks on the abscissa of Fig. 5 give the 

322 critical values x c for P{£ < x c } = 0.99, 0.999, 0.9999, and 0.99999 where P{£ < x c }is the probability 

323 that £ < x c . The critical value for P{£ < x c } = 0.999999 is just off the right of the chart at £ = 103.3 

324 Truncation of the pdf at any of the listed critical values will keep the integral in Eq. 30 finite and thus 

325 lead to a price for the option with the concomitant level of confidence. 

326 Truncation of the underlying pdf to obtain a price for an European option has raised some objections. 

327 It is felt that any price can be obtained with truncation and thus truncation is not a valid approach. 

328 Truncation is a valid approach. The standard model allows an infinite price, which is not physically 

329 possible. Truncation is one method to deal with a mis-specified model for the value of an asset. 

330 Another method is to use a model that takes into account a limited reservoir of money that is available 

331 to purchase an asset. The homogeneously saturated model is one such model. 
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332 Truncation also allows to quantify the risk inherent in setting a price for an option. If it is believed 

333 that the return on the asset follows a given pdf, then the writer of the option can select a probability 

334 that an asset will exceed a critical value and price based on the selected risk. This is no different than 

335 setting confidence intervals on experimentally determined numbers. Once the pdf for the experimentally 

336 determined numbers is known, the level of confidence is selected and this dictates the confidence interval. 

337 If one wishes for 100% confidence then the confidence interval will be infinite if the region of support 

338 is infinite for the pdf. If one is willing to accept some risk that the true value might lie outside the 

339 confidence interval, then the confidence interval is reduced. Conversely, one can select the confidence 

340 interval and determine the confidence level that is consistent with the selection. 

341 The logistic and homogeneously saturated models show indirectly that truncation (of the standard 

342 model) is a reasonable approach to determine the price of an option. These alternative models require 

343 knowledge of a saturation parameter f3 and the strength of the coupling between the money supply and 

344 the price of the asset. Knowledge of these parameters will be limited. Thus any option price determined 

345 from these alternative models will have uncertainty and a degree of arbitrariness associated with it. 

346 In addition, the noise driving term will need to be be limited or truncated to keep the money supply 

347 from approaching infinity. The uncertainty arising from truncation for the alternative models can be 

348 eliminated if the total amount of money available to chase the asset is known. However, this number 

349 will never be known with precision and will be somewhat arbitrary, much like the price of an option as 

350 determined from truncation of the pdf in the standard model. 

351 The homogeneously saturated model that is presented is this paper is not a sophisticated model but 

352 it does appear to predict returns that are consistent with returns that are observed. The model assumes 

353 no lower level population, as is appropriate for a simple III- V semiconductor diode laser [T7] . One could 

354 include a lower level, i.e, a reservoir of sellers, and transfer population or wealth between buyers and 

355 sellers. Depending on the relative sizes of the reservoirs and the strengths of the couplings between the 

356 reservoirs and assets, it might be able to mimmic operation of the markets and thus allow for realistic 

357 estimates of prices of assets on the market. 
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3 5 8 6 Conclusion 



359 In this paper, a homogeneously saturated model for the time development of the value of an asset 

360 has been considered and has been compared to the standard model and a logistic model for the time 

361 development of the value of an asset. These models were based on Langevin equations for the time 

362 development of the value of an asset. The logistic model uses a non-linear (quadratic) saturation term 

363 for the price of the asset. The quadratic saturation term is similar to the saturation term used in 

364 population models and in descriptions of Van der Pol's oscillator. The homogeneously saturation model 

365 borrows from lasers physics. Simple coupled rate equations were constructed and solved in steady 

366 state for the time development of the price about the steady state value of the asset and for the time 

367 development of the amount of money that is available to purchase the asset. 

368 The standard model is geometric Brownian motion and assumes normal statistics for the noise driving 

369 term in the Langevin equation that describes the time development of the price of the asset. It is known 

370 that returns have fat tails and that Student's ^-distributions describe well the returns for the DJIA and 

371 the S&P 500 indices. Hence simulations for the time development of the prices from the three models 

372 were performed with samples from a Student's t— distribution. 

373 Simulations showed that the logistic model limited the maximum return but not the minimum 

374 return. The homogeneously saturated model limited the magnitudes of both the minimum and maximum 

375 returns, and in this respect was, of the three models, the model that was consistent with observed returns 

376 for the DJIA and the S&P 500 indices. 

377 Both the logistic and the homogeneously saturated model require additional parameters, and in 

378 particular, a parameter that controls the (non-linear) saturation. The standard model does not allow for 

379 saturation and thus does not require a saturation parameter. It is likely that this saturation parameter 

380 will be somewhat arbitrary and uncertain; it will not likely be known exactly. Thus this lends some 

381 arbitrariness to the alternative models, similar to the arbitrariness that is involved in truncation of the 

382 underlying pdf in the standard model to obtain the price of an European call option when using the 

383 standard model with fat-tailed distributions. In this respect, the alternative models indirectly support 

384 the approach of pricing options by truncation of the underlying pdf in the standard model. 

385 The homogeneously saturated model, which is presented in this paper, is a simplified, one reservoir 
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386 model. The model could be extended to included multiple reservoirs to allow for the interplay between 

387 buying, selling, and parking money on the sidelines, and could be extended to include multiple assets in- 

388 teracting with the reservoirs. These interactions could be homogeneously broadened or inhomogeneously 

389 broadened, or a mixture of the limiting forms of the broadening mechanisms. 
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